NASA Technical Memorandum 105286 



Multiwavelength Pyrometry for Nongray 
Surfaces in the Presence of 
Interfering Radiation 


Daniel Ng 

Lewis Research Center 
Cleveland, Ohio 


-rs) 


. ■ .. . 


"IT V : T . v IT PY ; v y ; y 

" : T. ‘ ; •' ' 

e . . / 1 v CSCL 1 *♦ 


.5 ’! 



February 1992 




E-6617 


MULTI WAVELENGTH PYROMETRY FOR NONGRAY SURFACES 
IN THE PRESENCE OF INTERFERING RADIATION . * 

Daniel Ng 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

SUMMARY 

A NASA-developed multiwavelength pyrometry technique for nongray surfaces has been 
, extended to also measure surface temperature in the presence of interfering radiation. This 
^radiation is produced by heat lamps used to raise the temperature of the surface. The necessary 
instruments are a spectral radiometer, an auxiliary radiation source, and a computer. Four 
radiation spectra are recorded: (1) the unobstructed spectrum characterizing an auxiliary radi- 
ation source, (2) the unobstructed spectrum characterizing the interfering radiation, (3) the 
radiation spectrum consisting of surface emission plus the interfering radiation, and (4) a spec- 
trum consisting of the radiations of (3) plus the reflected radiation due to the incidence of the 
auxiliary radiation source on this surface. With these spectra, application of two-variable, 
nonlinear, least-squares, curve-fitting computer software determines the surface temperature and 
the spectral emissivity. Use of the method to measure the surface temperature of silicon carbide 
under a simulated interference condition is demonstrated at a low temperature just above 
ambient. The instrumentation necessary to extend the method to elevated temperatures is 
discussed. 

INTRODUCTION 

Pyrometry is a favorite noncontact surface-temperature measurement method. Tradi- 
tional pyrometers are the one- and two-color method. Multiwavelength pyrometers have, 
recently been introduced to measure surface temperature (ref. 1). , 

Ceramic materials are to be used in the development of the next generation of propulsion 
systems. Accurate surface temperature measurement using noncontact methods is urgently 
needed in the development of these materials. Most ceramics materials are nongray, and their 
emissivity is known to be wavelength and even temperature dependent. Emissivity is low at 
v short wavelengths, increases to almost unity at wavelengths around .8 or 10 ^m, and then varies 
, further at longer wavelengths (ref. 2). Therefore, traditional pyrometers are inherently 
inaccurate when applied measure ceramic surface temperatures, and unavoidable errors result 
(ref. 3). Furthermore, in the testing and development of these materials, intense external heat 
lamps are used as a heat source to raise the temperature of ceramic samples. The aforemen- 
tioned pyrometers cannot discriminate between the reflected radiations originating from the heat 
lamps and radiations emitted from the measured surface. 

Although we previously applied multiwavelength pyrometry to determine the surface tem- 
perature of gray surfaces in the presence of interfering radiation (refs. 4 and 5), that method also 
resulted in unavoidable error when applied to these nongray emissivity materials. 

Multiwavelength pyrometry was developed at NASA Lewis to measure the temperature of 
a nongray silicon carbide ceramic surface by simultaneously measuring the surface emissivity 


(refs. 6 and 7). The temperature measurement was accomplished with the acquisition and analy- 
sis of three radiation spectra from a 'nongray silicon carbide (SiC) surface. This multiwavelength 
pyrometry technique was thus extended to measure the temperature of a nongray surface in the 
presence of interfering radiation. 


EXPERIMENTAL METHOD 

Recall that\the basic components of this pyrometer are a spectral radiometer, an auxiliary 
radiation source which can be turned on and off (either manually or by mechanical chopping), 
and a computer for data acquisition, data processing and numerical calculation (refs. 6 and 7). 
The SiC surface from which measurements were made was heated by allowing it to equilibrate in 
front of a blackbody furnace. The SiC sample completely covered the blackbody cavity opening. 
The blackbody furnace temperature was regulated by a temperature controller with ±0.5 °C 
long-term stability. The auxiliary radiation source was regulated by a constant-current power 
supply such that its typical light output rms ripple was 0.05 percent. 

In addition to the three spectra required for the nongray multiwavelength pyrometry 
(ref. 6), an additional spectrum was also needed to characterize the radiation of the interfering 
radiation source. It was acquired by allowing the unobstructed radiation from the interfering 
source to enter directly into the spectral radiometer. The four required spectra were 

(1) S q (A): the spectrum of the auxiliary radiation source, which was acquired prior to real 
time measurement 

(2) Sj(A): the spectrum characterizing the spectral intensity of the extraneous interference 
radiation source 

(3) Sjj(A): the spectrum of surface-emitted radiation Sjj(A) plus interfering radiation 
Sjj(A) reflected from the extraneous source 

(4) Sjjj(A): the radiation spectrum Sjj plus the reflected radiation Sjy due to the 
incidence of the auxiliary radiation source on the surface 

Note that some of the definitions of the spectra S 0 , Sj, S jj and Sjjj are different from those 
designated by the same symbols in reference 6. 

The auxiliary radiation source was a 20-W infrared lamp. Spectrum S 0 (A) was obtained 
from the geometry shown in figure 1. S Q (A) needed to be measured only once for such a radia- 
tion source. 


Spectrum Sj(A) represented the spectral distribution of the interfering radiation. If this 
interference source were blackbody in origin, it would be represented analytically by a Planck 
function corresponding to its temperature. However, if it were more complicated in origin, its 
actual spectral distribution would be determined experimentally by using the geometry shown in 
figure 1. 


Spectrum S n (A) was obtained from a typical pyrometry measurement of a surface in the 
geometry shown in figure 2, with the auxiliary radiation source turned off. The emission spec- 
trum Sjj(A) of the nongray surface is described by a Planck function of temperature T modi- 
fied by the wavelength-dependent emissivity (ref. 6) plus the interference radiation Sjj(A) 
reflected by the measured surface. 
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Spectrum Sjjj(A) was also obtained from the same geometry shown in figure 2 but with 
the auxiliary radiation source turned on. Thus, Sjjj(A) was Sjj(A) plus reflected radiation Siv 
due to the auxiliary radiation source. Subtraction of spectrum Sjj(A) from Sjjj(A) produced the 
reflected spectrum Sjy(A) = Sj n (A) — Sjj(A) due to the incidence of the auxiliary radiation 
source on the measured surface. 


.'W By application of the arguments advanced previously (ref. 6), the measured reflectivity 
z(A) at any wavelength was a ratio given by 


z(A) 


Siv(A) _ S m (A) - S n (A) 



Because S fl (A) was measured on a direct path in figure 1 and Sjy(A) was measured on the 
reflected path in figure 2, the ratio z(A) included a constant factor f resulting from different 
conditions along the two paths. Two such conditions were the geometry of the optical beams 
and the nonspecular reflectivity surface. Thus, the true reflectivity r is 


r(A) = (2) 

f 


where f must be determined. 

From reference 6, the emissivity e and reflectivity r of a surface are e(A, /?, 0 , T e ) and 
r(A, /3, 6 , T e ), where and 9 are the polar angular coordinate angles specifying a direction 
with reference to a suitably chosen coordinate system and T e is the temperature of the reflect- 
ing surface. Even though the emissivity and reflectivity of a surface contain spectral, angular, 
and temperature dependence, only the spectral dependence remains important in pyrometry 
application (ref. 6). Therefore, by application of Kirchhoffs law, the emissivity of the surface is 
given simply by 

e(A) = 1 - r(A) r 

(?) 

f 

In principle, the constant f can be determined from the geometry of the experimental setup. 
However, it is more conveniently determined by curve fitting. 

The mathematical expression for the spectrum Sjj(A) is written as a sum of two terms: 

s„(A) - sJj(A) ♦ syj) (4) 

!'**> ^jj(A) describes the surface emission radiation. It is a Planck function modified by the 
rwavelength-dependent emissivity and is shown to be (assuming the direction of S 0 is not impor- 
tant or the same as Sj) - 
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exp 
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where c r and c 2 are the first and second radiation constants, and f and T are unknown 
parameters (ref. 6). 

Sjj(A) describes the reflected interfering radiation from the extraneous source Sj(A). It is 
proportional to the product of r(A), Sj(A), and a constant of proportionality g (due to beam 
geometry but different from f) 

S*,(A) = g r(A) 

Substituting equation (2) gives 

s n( A ) = g Z ( A ) 

If the ratio g/f is replaced by 

h = I 
f 

S^A) = h z(A) 

The quantity z(A) is obtained by using the wavelength-by-wavelength spectrum subtraction and 
division in equation (1). The adjustable parameters f, h (or g), and T are determined by 
applying-least squares curve fitting, with A and z(A) treated as independent variables to fit the 
spectrum Sjj(A) according to equations (4) to (7). 

The unobstructed spectrum S Q is shown in figure 3. It is the direct output of the detecor 
in volts. As described in reference 6, it can be easily converted into energy units using the 
response function of the detector. 

The interference source is simulated by a 40-W, nominal-temperature 644 K soldering iron 
which has a temperature controller of unknown temporal stability. It is placed in front of the 
SiC surface in approximately the same angular relationship as the auxiliary radiation source and 
is well shielded from the direct view of the spectral radiometer. 

The spectrum Sj of the radiation emitted by the interference source (in arbitrary energy 
units) is shown in figure 4. The other spectra necessary for temperature measurement are 
acquired using the geometry shown in figure 2. The surface whose temperature was to be meas- 
ured is a silicon carbide (SiC) wafer used previously (ref. 6). 


Si(A) 


Sj(A) 


Si(A) 


( 6 ) 


( 7 ) 
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. RESULTS . . . . . . 

The radiation spectra collected are shown in figure 5. Sjj(A) is the surface emission spec- 
trum at the test temperature. Sjj(A) is the emission spectrum obtained in the presence of the 
interfering source with the auxiliary radiation source turned off, and Sjjj(A) is the spectrum 
obtained with the radiation source turned on. From observation of Sjj(A), the contribution of 
the interference radiation to the pyrometer signal is clearly evident. At the short wavelengths, 
the interference component of the signal is over 50 percent. The total energy in Sjj obtained 
by integrating the spectrum from 2.5 to 14.5 /im is 1.13 times that due to surface emission 
obtained by a similar integration. If this contribution is not accounted for, the temperature 
determined will be in error. 

The spectra corresponding to Sjj and Sjjj in figure 5 are reproduced in volts in figure 6. 
The difference between these two curves is the reflection spectrum Sjy due to the incidence of 
the auxiliary radiation source S Q on the SiC surface (shown in fig. 7). Division of the spectrum 
of figure 7 by the unobstructed spectrum S Q of figure 3 gives the quantity z(A) defined by 
equation (1); the result is shown in figure 8. 

Commercially available nonlinear least-squares curve-fitting software was used to fit the 
spectrum Sjj in figure 5 to the mathematical expression of equation (4) with equations (5) 
to (7) explicitly substituted. Curve fitting yielded a surface temperature of 341.3 K. The values 
of f and h are 0.375 and 3.49, respectively. The fitted curve, the sum of the two terms in 
equation (4), is shown in figure 9 together with the measured surface emission spectrum. Also 
shown is the 341.3 K Planck curve. 

DISCUSSION 

To demonstrate the importance of correcting for the reflective radiation from, an extrane- 
ous source, spectrum Sj was omitted from consideration and only spectra S Q , Sjj, and Sjjj 
were subjected to the multiwavelength pyrometry analysis in which interference was not con- 
sidered (ref. 6). The temperature was determined to be 373.4 K. The 341.3 K temperature 
measured by the current multiwavelength pyrometer included correction for interference radia- 
tion but was not simultaneously verified by another method during the experiment. However, , 
from the good agreement between the fitted curve and the experimental spectrum shown in fig- 
ure 9, we assert that this is the correct temperature. It was therefore concluded that neglecting 
the effect of interference would introduce an error of about 10 percent. A 0.125-mm (5-mil) 
chromel-alumel thermocouple was subsequently attached to the sample to permit as much con- 
tact as possible between its junction and the silicon carbide surface. The surface was then 
allowed to equilibrate to a temperature under the same condition (including the presence of a 
simulated extraneous radiation source) as before. The thermocouple measured a temperature of 
339.7 K, which is in agreement with the multiwavelength pyrometry measurement to within 
1 percent. This result confirms the assertion that a good fit to the radiation spectrum is an 
indication of accurate temperature measurement. 

At higher temperatures, the surface emission will greatly overwhelm the radiation from 
the 20-W auxiliary radiation source. The reflected spectrum Sjy cannot be measurable by the 
subtraction of the two spectra obtained with and without the auxiliary source. However, the 
difference in signal can still be readily measured by chopping the auxiliary radiation source and 
using phase-sensitive detection schemes. Most spectral radiometers have an external signal 
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chopping capability. In this way, the subtraction of the two spectra, spectrum Sjj from spec- 
trum Sjjj, is automatically integrated into the signal processing and data acquisition as a single 
operation. The reflected spectrum Sjy is immediately obtained. After division by the direct 
spectrum S Q , straight forward application of two-variable, least-squares, curve-fitting computer 
software will determine the temperature of the nongray surface with the effect of interference 
automatically taken into account. 


CONCLUSION 

The multiwavelength pyrometry for nongray surfaces determined the temperature of a sili- 
con carbide surface in the presence of simulated interference radiation. Four radiation spectra 
were acquired and analyzed. The temperature was determined by application of two-variable, 
nonlinear, least-squares, curve-fitting software to fit the surface emission spectrum containing 
the interference radiation to a mathematical expression involving the wavelength-dependent 
emissivity of the surface and the interfering radiation source spectrum. The good agreement 
between an experimental spectrum and its fitted spectrum is an indication of the accurate 
temperature measurement that can be obtained by multiwavelength pyrometry. 
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Pyrometer Radiation 

source 

(auxiliary or interfering) 

Figure 1 . — Experimental arrangement to determine the unobstructed 
spectrum S 0 of the auxiliary radiation source or the interference 
spectrum S[. 
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Pyrometer Auxiliary 

radiation source 

Figure 2.— Experimental arrangement of pyrometry components. 
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Figure 3. — The unobstructed spectrum S 0 (in volts) of the auxiliary 
radiation source. 



Figure 5. — Spectrum of silicon carbide; S, 1 ,, surface emission; S„, 
surface emission plus interference (no auxiliary radiation); S n |, 
surface emission plus interference (with auxiliary radiation). 



2 4 6 8 10 12 14 

- Wavelength, pm 

Figure 4. — Spectrum Sj (in arbitrary energy units) of the simulated 
interference source. 



Figure 6.— Surface emission spectra (in volts) of silicon carbide 
surface; S[j(X) = sjft.) + S(j(X) (auxiliary source off); Sjjj 
( auxiliary source on). 
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Wavelength, pm 

Figure 7.— Reflected spectrum S[ V (In volts) from a silicon carbide 
surface obtained by subtracting Sjj from Sjjj. 



Rgure 8. — Measured reflectivity z{X) of a silicon carbide surface 
obtained by dividing Sjy by Sg. 



Wavelength, pm 

Figure 9. — Curves for spectrum Sij, spectrum according to 
equations (4) through (6), and Planck curve. Temperature, 
341 .3 K. 
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